Nonclassical rotational inertia for a supersolid under rotation 
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As proposed by Leggett (3J, the supersolidity of a crystal is characterized by the Non Classical 
f**"*. ■ Rotational Inertia (NCRI) property. Using a model of quantum crystal introduced by Josserand, 

Pomeau and Rica [5j], we prove that NCRI occurs. This is done by analyzing the ground state 
of the aforementioned model, which is related to a sphere packing problem, and then deriving a 
C^l . theoretical formula for the inertia momentum. We infer a lower estimate for the NCRI fraction, 

I ' which is a landmark of supersolidity. 
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PACS numbers: 67.80.-s 



A recent experiment by Kim and Chan 1] allowed them to measure the moment of inertia of solid helium and find 
that it is lower than its classical value. This property is referred to as Non Classical Rotational Inertia (NCRI). This 
experiment has raised a lot of interest and was interpreted as a landmark of supersolidity, on the basis of a paper by 
Leggett [|| . In Q , Leggett predicted that the property of nonclassical rotational inertia possessed by superfluid helium 
is shared by solids and proposed as a definition for the non classical rotational inertia fraction NCRIF = (Iq — I) / Iq 
where / is the moment of inertia of the crystal under study and Iq its classical value. One theoretical challenge 
(see the review paper of Prokof'ev Q) is to estimate this NCRIF and check that it is nonzero. This is the aim of 
this paper, based on a model of quantum crystal, introduced by Josserand, Pomeau, Rica [f|. In this respect, we 
derive a key estimate providing the lower bound © for the NCRIF. In the literature (see (H), different microscopic 
mechanisms have been proposed to describe the supersolidity of a crystal, based mainly on the off diagonal long range 
order property (ODLRO) of the density matrix, and Jastrow wave functions. Here, we do not relate directly the 
NCRI to ODLRO, or the presence of vacancies but choose another approach to model the solid. 

Josserand, Pomeau, Rica [f| proposed a model of quantum solid: it is based on the fact that the complex valued 
wave function common to all particles of mass m minimizes the Gross-Pitaevskii energy with an integral term that 
can be viewed as a 2-body potential in a first Born approximation: 

/ ^ |V ^ (r)|2 dr + \ jj ^ (r ' ~ r )\^( r ')\ 2 \iP( r )\ 2 dr dr ' 



where [/(•) is a potential depending on the distance between atoms. The normalization condition is v = J \ip\ 2 /V 
where V is the volume of the region T> occupied by the solid. This model bears an important difference with classical 
solids, in the sense that in classical solids, there is an integer number of atoms per unit cell, while in this quantum 
solid model, the average density is a free number, independent of the crystal parameters. Moreover, this model yields 
S* ■ a dispersion relation between the energy and momentum of elementary excitations that depends on the two-body 
potential. The choice of U is made in order to have a roton minimum in this dispersion relation. For instance, one 
possibility is to take f7(|r|) = Uo9(a — |r|), with 8(.) the Heaviside function ||. We define g = Uv^^-vc? and rescale 
distances by a so that the rescaled energy E g is given by 



||V^(r)| 2 d r + ^JJ U(r' - r)\^{r')\ 2 \ip{r)\ 2 dr dr' (1) 

where £/(|?"|) = 9(1 — \r\) and J \ip\ 2 = V. For small g, the ground state is \ip\ = 1, and for large g, computations in 
Q indicate the presence of a crystal phase with some supersolid-like behaviour under rotation. Moreover, the authors 
of H checked that this model also provides mechanical behaviours typical of solids under small stress. We believe 
that the model proposed in [5[ is not far from a realistic model of solid helium, that is of a dense solid with strongly 
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FIG. 1: The expected configuration of sets Ai in 2D. 

repulsive interaction. Note that for He, we have g ~ 25, and for Ne, g ~ 100 |6|. In the large g limit, we will see 
that the ground state of ([1]) is a periodic array of peaks. The self interaction of a peak becomes a constant, added to 
the energy, and independent of the wave function local profile. One could argue that in a mean field model of a real 
crystal, the interaction has a hard core, so that the self interaction is infinite. But in the true physical system of solid 
helium, a given atom does not interact with itself and thus does not provide any infinite self interaction. This, added 
to the various properties of the quantum crystal derived in [f| , which are in agreement with experimental solid helium, 
make us believe that the model provides insight into the understanding of supersolids. The aim of this paper is to use 
this model to derive an approximate theoretical value for the reduction of the moment of inertia of a supersolid. The 
proof is two-fold: on the one hand, we use the specific choice of the interaction potential U to get that for large g , 
the ground state ip g has a periodic density p g = tpg. Moreover, the wave function is localized around sets defined by 
a sphere packing problem. On the other hand, given this periodic density p g , we use the expression ([6]) of the NCRIF 
and the fact that tp g is a ground state, hence a solution of some nonlinear Schrodinger equation, to obtain a lower 
estimate of the NCRIF ([7J and ([5]) . Since this second part of the proof only relies on the periodicity of p g and the 
fact that it is a solution of an equation, it could be used for other models providing a periodic density, for instance 
that of an optical lattice (with an exterior field creating a periodic density with several atoms per site). 

If T> is a solid sample, the sphere packing problem 0,(1] provides a number n(V), the largest number of points in 
T> which are at distance larger than one from each other. When this number is large, the optimal location of points is 
proved 0,11] to be close to a hexagonal lattice in 2D. In 3D, 2 configurations are optimal: body centered cubic close 
packing and face centered cubic close packing. When g is large, the two terms in ([1]) are of different order, hence the 
ground state tp g is very close to a function iJjq that is found by minimizing the kinetic energy within the functions that 
minimize the interaction term, which is dominant. We are going to prove that such a function is supported in sets A4 
which are at distance at least one and whose number is n(T>). Thus, the sets Ai are determined by the minimization 
problem 




inf { > Xi(Ai) } (2) 

where Ai is the ground energy of —A in A4 with zero boundary conditions: Xi(Ai) = inf r |„|2 =1 { J A |Vw| 2 }. The 
expected configuration is illustrated in figure [TJ The function ipQ corresponds to the ground state of —A in each Ai 
and vanishes outside the Ai's. A ground state of E g will be very close to ipo in the sets A,, and exponentially small 
away from the A^s, except on a boundary layer. 

When the sample is set under rotation about the z axis, the free energy of the system is defined as 

e(n) = hrf{£ ff (vo-n<^,L,(vO)} (3) 

where L z (ip) — irxWip-e z and E g is the energy defined in |T]). When is small, e(Q) can be expanded as eo — (l/2)/f2 2 
where I is the effective moment of inertia of the system. Leggett [4] suggested as a criterion for superfluidity the 
existence of a non classical rotational inertia fraction (NCRIF) , defined as (Jo — J)/ Jo ; where Iq is the classical moment 
of inertia of the crystal phase and is equal to J |^ 9 | 2 r 2 where ipgis a ground state of E g . The point of this analysis 
is to find an estimate for the NCRIF, computed numerically in [5J, and prove that it is non zero. As can be seen in 
([6]), a good knowledge of p g , the ground state density with no rotation, is needed to estimate NCRIF in ([7]) and |(5J). 

The paper is organized as follows: first, we study the ground state of the crystal phase with no rotation and derive 
@. Then, we present some more refined computations in the ID case, and finally we derive estimates for the NCRIF. 
Crystal phase with no rotation We first describe the minimization of the second term of |T]) which provides a class 
of functions ip such that p — \t[)\ 2 has mass located in disjoint sets Ai, at distance at least the range of the potential, 
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FIG. 2: Ground state of E g when g is large (dashed line), and its limiting profile ipo (solid line). The bumps are of size I and 
separated by a distance I. 

which is 1. When one wants to minimize / | Vip\ 2 in this class, this provides a constraint @ on the shape of the sets 
Ai that we explain. 

We denote by (U * p)(r) = J U{r - r')p(r')dr' and F{p) = J J U{r' - r)p(r')p{r) dr dr' . Recall that n(V) was 
defined in the introduction and is related to the sphere packing problem. 

Theorem 1 A measure p with J p — V minimizes F(p) if and only if there exist n(T>) pairwise disjoint sets 
A%, . . . , A n (p-s, such that 

f V 

dist(Ai,Aj) > 1 ifi^j, and / p= ——. (4) 

Ja z n{T>) 

Moreover, mini* 1 = V 2 /n(T>). 

The proof of this result which strongly relies on the shape of U is made in the appendix. Let us call po a ground 
state of F and ip g a ground state of E g with p g = ip 2 . Then F(p ) < F(p g ) < F(p ) + (l/g) J \Vip \ 2 - |V-0 9 | 2 . For 
g large, we deduce that p g is an almost ground state of F. Among all the possible po's which are ground states of 
F, the limit of tp g when g is large should be such that ipo = ^fp~o minimizes the kinetic energy J |V - 0| 2 among all ip 
such that p = |^| 2 is a ground state of F. This implies that the support of ipo is the union of n connected sets Ai, 
1 < i < n, which satisfy (J5J). 

More specific computations in ID. In dimension 1 (for N atoms in a cylindrical annulus Q), that is if V — (0, L), 
then n = n(V) = [L] + 1 (if L is not an integer [l(|), and the Ai's are intervals (xi,Xi + I), with I = L/n — 1 + 1/n 
and Xi — i(l + 1). Thus, ipo{x) = (2L/nl) sin(7r(a; — Xi)/l) if x S (xi,x i+ i) and otherwise (see figure[2l). Moreover, 
Eo(ipo) = tt 2 L/21 2 . Indeed, the ground state of F provides n sets Ai separated from one another by distance at least 
1. Hence Ai C [aj,6j] for all i, and hi + 1 < aj+i for i = 1, . . . , n — 1. Then, since ([J} implies that f ' \u\ 2 — ~ for 

every i, J \u'\ 2 = \u'\ 2 > ^Li la, \ u \ 2 = ^ E"=i ( bi - ai )^ with quality if and only if the restriction 

of it to each interval (dj, bi) is a scaled and normalized sine function multiplied by a constant of modulus 1. Moreover, 
Jensen's inequality implies that ^2 npj^-p > ( n -i ^Jj>--a-)) 2 — n /^ 2 ' equality if and only if b t — a, = I for every 
i. 

One expects a boundary layer around each Ai . In this one-dimensional setting, it is possible to compute it explicitly. 
In order to do so, we assume that ip is a dilation of the limit ipo, namely ip(x) = (i+k)n sm ("^"T+F"^) ^ 
x 6 + 1) — k/2, i(l + 1) + 1 + k/2], and ip = otherwise. The energy of this trial function is computed explicitly, in 
the limit of small k: E g (ip) w + + jA (j^) , where A = 13L 2 7r 6 /(90n). Minimizing this expression with 

respect to k yields k = ((2L7r 2 Z 3 )/(3A)) g^ 1 ^ 5 . Inserting this into the expression of the energy, we find 

„,,s * 2 L gL 2 5 /60L 4 7T 6 n 2 \ 1/5 _ 1/5 

as g is large. 

The above computation indicates that the boundary layer around each bump of the limit function ipo is of order 
g" 1 / 5 and if x denotes the scaled distance to the boundary, then the matching between ipo and in the boundary 
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layer is described by the solution of u" = cx 3 u. This boundary layer decreases the energy by an amount of order 

Dimension 2 and 3. There is no complete determination of the A^s, except that once the sphere packing problem 
is known to provide a hexagonal lattice, the Aj's are sets whose centers are located on an almost hexagonal lattice. 
Since minimizing Ai(A^) over Ai with fixed volume implies that Ai is a ball (see Q), condition ([2]) implies that Ai 
"looks like" a ball to some extent. However, Ai(Aj) is increasing with respect to Ai, which implies that Ai cannot be 
exactly a ball, but is closer to a hexagon (see figure [T]). 

A ground state ipg of E g is close to i[iq in the Ai 's and small in between. We need to understand better the smallness 
of Vv The Euler-Lagrange equation satisfied is — Aip g + gW(x)ip g = 0, where W(x) = U * \ip g \ 2 — X g , and A g is the 
chemical potential. The shape of U and the mass constraint on tjj g imply that 1 < maxIW^I < 2. Thus the Harnack 
inequality applied to the equation for tjj g yields 

inf ip g > 2gC d e ( - V ^ T) max^, (5) 

for some constant Cd [HJ , where T is the size of the period of p g and is of order 1. This estimate is used below in 
order to estimate the NCRIF. 

In the limit of very large g, the function ip g is exponentially small: the potential W is almost equal to Wq — U *po— Ao, 
which vanishes in each Ai, and is positive in between. Using appropriate comparison solutions, it is possible to prove 
the estimate ^(cc)! < exp(— S^/gas), for any x such that dist(a;, Ai) > 8. The constant as is the minimal value of W 
in this region, and is of order S^ d+5 ^ 2 , where d is the dimension. The density is thus exponentially small between the 
Aj's. However, in the experiments, g ~ 25, so that it is not large to the point of having tiny density. 
Small rotation When the sample is set under rotation SI about the z axis, the free energy of the system is defined by 
^ and E g is the energy defined in (TTJ). We assume that the ground state ip of §5§ is of the form ip(x) = p g (x)e tns ^ x ' 
for small SI, where ip g = ^fpg' is a ground state of E g , that is for SI = 0. This equivalent to expanding the phase in 
terms of SI and assume that the first order variation in the phase is not sensitive to the variations in density in terms 
of SI. Then, the phase S should minimize J p g \V S — e z x r| 2 among all possible test functions. This provides an 
expansion of e(Sl) for small SI and hence a value for / which allows to compute 

inf a f PolVS* — e z x r| 2 
NCRIF = — bJ F9 \ s—? L (6) 

J Pg r 

Two limiting cases are easily identifiable: when p g = 1 (i.e when g is small), this ratio is 1, and when p g is periodic 
and has all its mass localized in the center of the cell, this ratio tends to 0. For intermediate values of g, the wave 
function is localized in Ai, with tails in between the sets which are small, but not too small. Then, ([6]) implies 

NCRIF > mf ^ ^/|Vg-e,xr|» = mfp^ 
max p g J r 2 maxp s 

The last equality is due to the fact that infg J \VS — e z x r| 2 is achieved for V5 = 0. Note that this estimate is a 
mere consequence of ([5]) and not of the shape of p g . The ratio maxp g / inf p g was estimated above ([S]). We thus have 

NCRIF > 4g 2 Cje~ 2V2 ~° T , (8) 

Since ^fg is of order 5, this implies that NCRIF ^ for the experimental values. 

In the very large g limit, ([5]) is not so good, but in this case, we may replace p g by po in Moreover, in each 
Ai, we can define local coordinates rj with respect to a point in Ai whose coordinate is x«. Then r = + x« and the 
phase S can be defined as a local phase Si in each Ai through VS = VS* + x^" where x^=e 2 xxj. We thus have 

NCRIF » [ni s ' ^ Po ^ Sl ~ GzX 1-1,2 



J Por 2 

Assuming that each Ai is the translation of a reference set Aq, the numerator is proportional to n(T>) times the infimum 
of the cell problem, which is always less than V VoI(Aq). Note that this cell problem depends on the volume since the 
size of Ao depends on n(T>). If V is large, a coarse-grained approximation for po yields that j p$r 2 w J r 2 cx V^ d+2 ^ d , 
where d is the dimension bigger than 2. Hence NCRIF < V~ d/2 Vol(A Q ), which tends to in the limit of large 
V. However, according to Legget [H, the system can be considered as superfmid if NCRIF ^> 1/N, where N is the 
number of particles equal to vV) where v is the initial average density (included in our rescaling providing g) . In a 
thermodynamic limit with V large g large, and v not fixed but large as well, we may still have NCRIF ^> 1 /N. 
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Let us point out that this behaviour contrasts to the ID case, where, in the largeo asymptotic, the NCRIF is zero: 
indeed a similar computation yields that it is equal to L 2 /(J p g J l/p g ) (see also [J]). In the large g case, this tends 
to since p g tends to po which is compactly supported and thus J 1/po = oo. 

Conclusion: We have derived properties on the density of the ground state of a quantum crystal. This has allowed 
us to estimate the NCRIF and find that the system is supersolid, on the basis of a definition of Leggett [4|. This 
complements the results of [5| and provides theoretical justification of the non zero NCRIF in the experiments of Kim 
and Chan 

Appendix: Proof of Theorem 1. Step 1. If p satisfies ((4]), we prove that F(p) = V 2 /n. Indeed, for every j, 
if x & Aj, then B x n (UjJ =1 Ajt) = Aj, where B x is the ball of radius 1 centered at x. Since Q implies that 
p(T> \ LiAj) = 0, we get that if x € Aj, then p(B x ) = p(Aj) = V/n. Since UAj is a set of full p measure, we get that 
F(p) = J(U *p)(r)p(r)dr = Ep 2 (^) = V 2 /n. 

Step 2. Let p be a ground state for F. We can argue by induction that there exist n = n(T>) points x%, . . . ,x n such 
that 

\x. t - Xj\ > 1 and U* p{x t ) = inf U * p = F(p)/V. (9) 

The last equality is in fact the Eulcr-Lagrange equation. The definition of n implies that if the Xi are any such points, 
then L)B Xi D V. So that V < p(U t B x J < £\ p{B Xi ) =Y ll u * p( x i) = V F (p)- Tllus m[nF = F (p) = y2 / n - 

Step 3. We have to check that ([3]) holds. For each Xj, we define Aj = {x £ B Xj n suppp}. Then by the Eulcr- 
Lagrange equation, U * p(x) — V/n in Aj, hence p(Aj) — V/n. We have that /j(UAj) = p(Ai), hence p(Ai PI Aj) = 
and if yi £ Ai, then U * p{yi) — p(B y . ) = V/n. Hence the points ({xi, . . . , x n } U {yi}) \ {xi} satisfy ^ and this proves 
©• 
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